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PREFACE 

This  report  deals  vlth  tho  necessarv  first  nhase  of  analvsls  of  weather 
T»dlflcatlon  experlnents,  namely  the  process  of  summarization  of  precipitation 
data.  Msually  precipitation  attrlhutahle  to  experimental  units  Is  measured 
through  a network  of  ralnyap,es.  'lethods  of  summarlzlnf;  such  data  can  vary 
In  sophistication  from  the  use  of  simple  averages  to  asgr'^^^ated  responses 
obtained  from  a response  surface  or  throuph  use  of  ma1or  components  of  vari- 
ation (say,  principal  components). 

In  this  report  v;e  examine  the  use  of  response  surface  flttlnp,  in  appli- 
cation to  data  from  Phase  I of  the  Santa  Barbara  Convective  Band  Seeding 
Test  Program  conducted  for  the  'law  hy  North  American  Heather  Consultants. 

"e  are  Indebted  to  NAHC  for  provision  of  data  tapes  for  our  analyses. 

A second  technical  report  on  the  use  of  principal  components  In  pre- 
cipitation data  summarization  Is  anticipated.  Other  technical  reports 
associated  with  this  contract  are  listed  at  the  end  of  this  report. 


Ralph  A.  Bradley 
Principal  Investigator 


T>AT.^  rt^^^'AFIZATTO"  T’!  A ’'nnil^TCATiniJ  Fy’’'’UT**'7’T; 

I.  A ’?.FSPn”Sr  ‘',Tiri?AC»'  APPPOACH^ 

hv 

Ualph  A.  PraHley,  SusMl  S.  Srlvastava  and  Adolf  Lanzdorf 
Department  of  Statistics,  Florida  State  Unlversltv,  Tallahassee,  Florida 

1 . INTRODDSTTO^’ 

The  fundamental  and  first  question  in  rep.ard  to  demonstration  of  seedlne, 
effects  Ip  weather  modification  experiments  to  enhance  precipitation  relates 
to  appropriate  summarization  of  precipitation  data  from  an  experimental  unit 
of  such  an  experiment.  'Experimental  units  In  weather  mo<llf icatlon  experimen- 
tations are  usually  either  a fixed  unit  of  time,  a alven  number  of  hours,  a 
dav,  a season,  or  alternatively  the  duration  of  some  storm  system.  Variation 
over  experimental  units  mav  he  so  large  as  to  forestall  any  mraningful  study 
of  treatment  (seeding)  effects. 

One  of  the  objectives  of  the  research  on  the  design  and  analysis  of 
'iTcather  modification  experiments  at  the  Florida  State  'Jnlverslty  is  the 
development  of  better  statistical  methodologies  to  summarlza  and  measure 
precipitation  data  from  an  experimental  unit.  If  this  can  he  done,  a 
contribution  ’'111  he  made  towards  better  formulation  of  statistical  models 
and  Improved  precision  In  consideration  o'  treatment  (seeding)  effects. 

Notwithstanding  the  Imnortance  of  summary  measures  of  precipitation, 
the  weather  modification  literature  contains  little  on  efforts  In  this 
direction.  Host  studies  use  simple  averages  of  ralngage  measurements  or 
at  least  some  function  of  these  averages,  for  example,  ratios  of  target- 
to-control  Area  means.  Since  ralngage  measurements  represent  rainfall  at  a 
set  of  points, each  of  which  represents  a much  larger  and  Ill-defined  area. 
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averages  based  on  thein  or  any  transformations  thereof  may  not  give  a good 
renresentatlon  of  true  total  rainfall  In  anv  desired  area.  There  Js  clearly 
a need  for  Investigation  of  improved  summarization  methods,  oarticularly  if 
the  objective  is  to  Increase  precipitation  in  the  ^ratershed  of  a reservoir. 

Research  at  the  Florida  State  University  has  the  objective  of  the 
development  of  statistical  methodology  for  Improved  design  and  analysis  of 
weather  modification  experiments.  But  applications  of  new  methodology  ‘ite 
necessary  and  the  NAUC  3anta  Barbara  data  was  made  available  for  this  purpose. 
Initial  work  under  this  contract  was  divided  into  two  main  parts,  summarization 
of  precipitation  data  for  each  experimental  unit  and  the  Identification, 
summarization,  and.  In  some  cases,  recovery  of  concomitant  information  on  the 
experimental  unit  that  mav  correlate  vrlth  precipitation  and  reduce  the  Inherently 
high  experimental  errors  associated  with  cloud-seeding  experimentation.  Our 
basic  t^net  is  that  Improved  experimental  design  can  best  be  obtained  through 
Improved  Identification  and  measurement  of  concomitant,  cloud-physics  variables. 
This  report  deals  with  surmarlzation  of  precipitation  data.  Other  reports  will 
deal  with  available  concomitant  variables  on  the  Santa  Barbara  data.  Use 

of  covariance  analvsls  for  reexamination  of  the  effects  of  cloud  seeding  is 
planned  when  appropriate  data  summarization  is  complete. 

In  this  report,  major  attention  is  focused  on  the  use  of  response  surface 
methods  to  summarize  precipitation  data.  Other  approaches  Included  the  use 
of  multivariate  methods  and  examination  of  ratios  of  average  target- to-control 
area  precipitations.  Some  revlevr  of  the  Santa  Barbara  experiment  is  given. 

This  report  is  broadly  organized  in  the  seouence  in  i^hlch  the  work  proceeded. 
Since  the  method  of  sunr.iarizatlon  proposed  is  tested  with  data  made  available 
by  NAWC  from  Phase  I of  their  Santa  Barbara  Convective  Band  Seeding  Teat  Program, 
Section  II  deals  with  a brief  description  of.  these  experiments  and  the  data 
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penerated  from  them.  The  analyses  of  the  generated  data  as  done  by  NAWC  and 
their  conclusions  are  also  contained  In  this  section,  Section  ITI  presents 
a discussion  of  data  used  In  response  surface  study  and  their  comparisons  and 
contrasts  with  data  used  in  two  reports  prepared  by  HA'TC.  A pilot  study  to 
determine  the  appropriate  functional  form  of  response  surfaces  was  first  Initiated 
through  use  of  observations  from  twelve  experimental  units.  A discussion  on 
the  alternative  forms  of  response  surfaces  fitted,  analyses  of  distributional 
patterns  of  residuals  from  these  surfaces,  and  the  distrlbutlr nal  patterns  of 
the  observed  precipitation  data  are  contained  In  Section  TV.  Section  IV  also 
includes  discussion  of  a basic  computational  problem  due  to  badly  centered 
independent  variables  in  regression  models  that  arose  In  the  Initial  oh.ase  of 
the  surface  fitting  work.  The  pilot  study  showed  that  the  appropriate  forms 
of  response  surfaces  were  two-dimensional  cubic  functions  with  location 
coordinates  as  Independent  variables.  The  cubic  surfaces  were  fitted  to 
ralngage  measurements  for  control  and  target  areas  separately  and  for  each 
experimental  unit.  A discussion  of  these  surfaces  as  representative  of 
precipitation  response  data  is  contained  In  Section  V.  Finally,  Section  V 
contains  the  '’er<vat5cp  cf  summary  Indices  of  precipitation,  representative 
of  rainfall  volumes  in  designated  areas  and  their  relationships  v^ith 
simple  averages  of  ralngage  measurements. 
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IT.  T)IE  SANTA  BARBARA  RXPRPTrRT'’T  A^’P  DATA  ANALYSTS 


2.1  The  Background . '^allent  features  of  the  Santa  Barbara  Convective 
Band  Seeding  Test  experiments  and  their  data  analyses  as  reported  bv 
Elliott  and  Tliompson  and  Thomoson,  Brow,  and  Flllott  are 
suTTnarlzed.  The  summary  Is  Intended  to  give  an  overview-’  of  the  ex- 
periment and  analyses  to  provide  a meaningful  background  to  the  data  set 
used  In  our  suggested  precipitation  data  summarization. 

The  Santa  Barbara  experiments  ’-'ere  conducted  In  t'^'O  phases  between 
'he  years  1°A7  and  1^74  by  North  American  Weather  Consultants  under  the 
''nsorship  of  the  ^'aval  '.Teapons  Center  (*”'C),  China  Labe,  California. 

Phase  I of  the  experiments  was  conducted  from  the  1Q^7-AR  season  to 
the  197A-71  season';  Phase  IT  was  from  1971-7?  to  1^73-74.  These  ex- 
periments had  the  oblectlve  of  demonstrating  the  effectiveness  of  cloud 
seeding  and  seeding  methods  In  west  coast  cvclonlc  winter  and  spring 
storna>  During  Phase  I,  the  primary  seeding  mode  »Tas  ground  based  (from 
a mountain  crest  at  1065  meters  above  mean  sea  level)  with  use  of  a 
hlgh-output,  silver-iodide,  pyrotechnic  device  Ignited  at  intervals 
during  the  passage  of  convective  bands  (defined  later),  ’'urlng  P^ase 
II,  the  seeding  mode  was  aerial  with  a continuously  burning  acetone  - 
Ag  I - NH^I  jet  seeder.  Bone  experimental  air-seeding  was  done  concurrentl'V 
with  Phase  I oneratlons’  the  results  were  Included  In  Phase  TT  analyses. 

Some  ground  based  seeding  was  done  during  Phase  TT  also,  mainly  as  a 
back-up  for  aerial  seeding.  The  focus  In  this  section  is  on  precipitation 
data  from  Phase  I since  these  data  are  used  in  application  of  the  response 
surface  methodology.  Later  analyses  of  Phase  II  data  are  planned. 

All  of  the  raingage  stations  which  were  upwind  or  west  of  the  ground 
seeding  station  were  designated  as  control -area  stations  and  the  ones 
downwind  or  east  of  the  seeding  site  were  designated  as  target-area  atatlons. 
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The  exnerlmental  unit  for  the  Phase  T "^anta  Barbara  exnerlments  was 


chosen  to  he  a convective  han^l  - a series  of  clouH  cells  arranged  In  a 
line.  This  choice  "as  node  on  the  basis  of  Investigations  by  Aeronetrlc 
Research  Tnc.  (an  affiliate  of  'TAT'o)  under  contract  to  the  ’’ational 
Science  Foundation.  ’Procedure  for  Identification  and  tracMne  of  convective 
hands  ^ras  developed.  This  "’"ocedure  consisted  of  alerting  the  seeding 
technician  of  a probable  hand  passage  If  either  (ll  one  of  six  stations^ 

Ir.  the  control  area  registered  at  least  '1.02"  rain  In  a fifteen 
minute  nerlod  or  (11)  the  radar  operated  by  reported  a banded  echo. 

Confimation  of  band  passage  was  given  In  case  of  an  alert  through  pro- 
cedure (1)  above  If  either  another  of  the  six  control-area  stations  registered 
a 0.02"  rain  In  a subsequent  fifteen  minute  nerlod  ...  or  the  radar  operator 
reported  an  incoming  convective  hand.  In  case  of  'alert'  through  procedure 
(11),  the  band  was  confirmed  If  any  of  the  six  control-area  stations 
registered  0.02'  rain  In  a fifteen  minute  nerlod  provided  this  coincided 
with  the  radar  position  of  the  band.  Further,  the  seedabllity  of  a con- 
firmed banfl  was  '’etermined  on  the  basis  of  wind  flow  compatible  with  the 
expectation  that  etfects  of  seeding,  would  fall  mainly  in  the  target  area. 

A.n  additional  seeding  criterion  based  on  air-mass  structure,  was  specified 
but  could  not  he  used  In  actual  operation. 


^These  stations  are  '’3,  ‘^^Ol,  ‘5206,  S235,  5:251, 


and  0 7^40. 


Table  1 below  j’.ives  the  distribution  of  bands  by  seasons  and  by  months 
within  each  season  for  the  four  seasons  of  Phase  I operations.  It  Is  seen 
from  this  table  that  In  all  107  bands  were  considered  seedable.  Actually, 
several  more  bands  considered  seedable  during  January  and  February  of  the 
19f’J-69  season  were  not  seeded  since  these  months  turned  out  to  be  months 
of  unusually  heavy  rainfall  and  seeding  operations  were  terminated 


Table  1:  DISTRIBUTION  OF  BY  SEASONS  AjTD  f'CNTTr.  : 

P’JASF  T 


Season 

f’onth 

NOV. 

PFC. 

JAN. 

FEB. 

’'lAR. 

APR. 

TOTAL 

1967-AP. 

i 

- 

f(7) 

8(3) 

4(3) 

1 

1 

4(3) 

2.3(11) 

1963-69 

4(1) 

6(?) 

17(9) 

11(5) 

UO) 

2(1) 

41(1'') 

1069-70 

0(1) 

- 

c.(A) 

5(4) 

6(4) 

- 

22(13) 

1970-71 

7(5) 

10(5) 

- 

2(1) 

1(1) 

2(1) 

22  0 3) 

13(7) 

16(8) 

32(15) 

26(13) 

12(3) 

3(5) 

107(56) 

Flpures  in  parenthesis  are  numbers  of  bands  seeded. 

Source:  Constructed  from  '7  , Appendix  C. 

dvie  to  threats  of  flooding,  A random  process  was  used  to  select  56  bands 
for  seeding.  The  project  was  Initiated  In  mid-January  of  1963  and  therefore 
there  are  no  data  In  the  table  for  the  months  of  November  and  December  of 
the  first  season.  The  largest  number  of  seedable  bands  occurred  In 
January  and  February  of  1968-69 • these  were  months  of  above-average  monthly 
rainfall.  March  and  April  for  all  the  four  seasons  and  November  and  December 
for  1969-70  had  below-average  monthly  precipitations.  Further,  there  vras 
a drought  in  Southern  California  from  January  In  1971  and  thus  there  were 
only  5 seedable  bands  In  as  many  as  four  months  of  the  1970-71  season. 
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2.2  Ha.qlc  Information  Collected  nnd  '’road  Out  lines  of  Analvses  ”nde . Durinr^ 


the  passage  of  a band,  a nunber  of  measurenents  were  recorded.  The  Santa 
Barbara  data  are  amonp  the  best  available  since  attention  was  given  to 
statistical  design  and  a substantial  amount  of  concomitant  information  was 
collected.  Important  information  gat  .'red  during  r!ie  course  of  band  passage 
is  as  follows' 

(1)  ilaingage  perclpi  tat  ion  measurements  from  a netvrork  of  stations  in 
control  and  target  areas. 

(il)  Estimated  duration  of  band  passage  ume  in  minutes  at  each  ralngage 
station . 

(Hi)  ;Junber  of  fusees  ignited  by  band. 

(Iv)  Alrnass  structure  of  tbe  band  in  terms  of  one  of  three  stahilltv 
categories . 

(v)  500  mb.  temperature. 

The  network  of  stations  for  which  the  estimates  of  total  precipitations 
and  durations  hy  hand  wore  obtained  is  sl'iilur  to  that  of  th<>  -ap  nr^- 
sentod  in  'pponiily  ''1,  ure  '.-I.  Airmass  structures  and  500  mb.  temperatures 
were  obtained  from  radiosonde  soundings  with  G!'n-1  equipment  at  S.anta 
Barbara  airport  along  wltb  similar  soundings  taken  at  nearby  government 
installations  during  the  course  of  band  nassafe.  Since  'start'  and  'end' 
times  of  band  passage  rrere  noted  ' s ' and  ' 7 , the  Information  on  whether  the 
band  passed  through  in  night  or  dav  could  also  he  derived. 

Precipitation  and  duration  analvses  were  ilo;ie  using  what  were  termed 
composite  single  ratios  (S.1.)  and  composite  double  ratios  (D.R.).  The 
•word  'composite'  is  used  to  indicate  that  these  ratios  use  data  from  all  the 
four  seasons  in  contrast  to  Individual  season  data.  Calculation  procedures 
for  S.R.  and  P.H.  are  explained  below  in  terms  of  precipitation  data. 
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L t denote  nr  ic  initation  at  fitatim  i f»"nri  band  bet  = 1 nr 

ij  ■ -i 

as  band  j Is  seeded  or  not  seeded  and  station  1 is  operative,  6^(i)  = 

the  nunber  of  seeded  bands,  while  H(i)  - M ''i')  * (t)  'l-^.Cl)  i.  the 

s ' ns  ' ^ .T 

nunber  of  non-seeded  bands  recnr  'ed  at  station  1 , (1)  b?ln;;  the  total  number 

of  seoichlv  ‘•■-.-c's  r_  -ordvd.  Precipitation  averap.es  were  defined  as 


f = y a ''l''v  /N  (1) 
s(l)  J Ij'  o' 


and 


ns  (i) 


These  are  station  averap.es  for  station  1 respectively  for  seeded  and 
non-saeded  '^ands.  bet  T index  the  six  control-area  stations  used  In  band 
detection.  Then 


b«>  ■ J 


and 


bs«)  ■ f 


and  C_  and  are  the  corresoondinp  averapes  over  the  available  control 
area  detection  stations.  For  station  1,  the  composite  single  ratio  was 
defined  as 


F.R.(l)  - T ,.JT  ... 

s(l)  ns(l) 


(1) 


and  the  composite  double  ratio  was 


n.R.d)  - (f  ,.v/0-)/(f  ,,v/C--), 

s(l)  s ns(t)  ns 


(?) 
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Tt  is  stated  that  each  of  these  sets  of  composited  single  and  double 
ratios  was  subjected  to  a * 'ann-'^il tnev  U test  to  test  the  null  hypothesis 
that  seeded  and  un-seeded  precipitations  were  dravm  from  the  sane  population 
versus  the  alternative  that  precipitations  from  seede<*  hands  were  stochastically 
larger  than  those  from  the  unseeded  hands,  "ann-’.'hltnev  tests  were  nerformed 
for  each  station  separately.  Definition  of  the  relevant  tv/o  samples  Is  not 
clear  from  4,  7'.  It  annears  that  the  test  of  significance  at  station  1 
for  the  single  ratio  is  simply  based  on  the  precipitation  measurements  v^^ 
separated  into  tv;o  sam'iles  as  band  1 Is  seeded  or  not.  In  regard  to  use 
of  the  composite  double  ratios,  the  precl pltations  v^^  may  have  been  used 
in  ratios  y^j/C^  v'here  Is  the  precipitation  average  for  the  (at  most) 
six  control  stations  for  band  1 , the  division  of  these  ratios  Into  two 
samples  being  again  on  the  basis  of  seeding  or  non-see/ ’ •’g . It  Is  noted 
that  comnosite  double  ratios  cannot  ho  obtained  from  Individual  members 
of  the  two  samples  as  defined  above. 

Double  ratios  tfere  designed  to  Include  correction  fur  the  natural 
Intensity  of  a given  banct  throug.h  division  by  a measure  of  control-area 
average  precipitation.  However,  various  analyses  indicated  that  normalizing 
the  single  ratios  (1)  hv  the  use  of  control -area  averages  to  obtain  double 
ratios  did  not  materially  change  the  results.  Therefore,  it  war.  determined 
that  a control  was  not  needed  for  evaluation.  Thus,  the  final  report  '’7'' 
contains  analyses  only  of  single  ratios. 

Table  2 gives  t^ie  distribution  of  stations  by  the  magnitude  of  single 
ratios  for  precipitations  In  Phase  I of  the  .Hanta  Barbara  experiments. 

This  table  also  contains  the  breakdown  of  this  distribution  betv'een  target- 


and  control -area  stations. 
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Table  2-  niSTRI JUTIO*!  0^  NIIMBFR  ttaTTOMR  bY  T?4jinrvS 
or  IIMCLT  TATIOS  T PB'TTPTTATIO’T 


RANGE 

i.P95 

and  loss 

1 .'>^0  - 1 
1.1P5 

.203  - 
1.3Q5 

1 .400  - 
1 .595 

l.bOO  - 1.300  - 

1.705  I.005 

2.O00+ 

I”  T0Ty>L 

Of 

TAPGET 

1 

e 

20 

22 

7 

1 

4 

01 

STauoNS 

CONTROL 

6 

20 

10 

O 

- 

- 

30 

ALL 

7 

2b 

-•-r-ii-  -T 

30 

25 

7 

1 

4 

Ipo 

Source.  I)erived~rron.  'X7  ^pendix  TI].  pnpe  ^ for  definitions  of  Tarpct  and 

Control  see  Table  'i  bolo"  for  d<»fl.nltlonn  unod  in  our  dot"'  sjirmnriznrlnn 


It  Is  seen  from  this  Table  2 that,  ''bile  t*^e  mean  of  slnp.lo  ratios  Tor 
target-area  stations  anpears  h^'’her  than  that  for  control-area  stations, 
both  target-and  control  area  averages  exceed  one.  '’’his  Is  Indicative  oT  higher 
precipitations  from  seeded  bands  than  non-seeded  bands.  But  it  happens  also 
in  control  areas,  a posslbllltv  suggesting  contamination  effects.  There  are 
27  out  hi  target  area  stations  which  recorded  values  of  F.R.  less  than  1.4. 
''ulte  a fe’'  of  the  22  stations  In  tlie  next  '’.P.  class  of  1.4'’'n  to  1.59S  had 

« 

In  frct,values  less  than  l.S.  Therefore,  about  half  o'"  the  target-area 
stations  mav  have  had  precipitation  Increases  from  seeded  hands  in  excess  of 
50  percent.  It  Is  noted,  though,  that  4 stations  in  the  B.P.  class  .2.00'Vf 
seem  to  be  outlying  observations.  Blngle  ratios  In  then  are  unduly  large  - 
one  as  large  as  7.32  (Station  A27),  while  a neighbouring  station  A23  has  the 
smallest  observed  S.R.  value  of  onlv  '>.79. 

Caution  in  Interpretation  of  Table  2 is  reoulred.  'Piile  effects  of 
seeding  seem  indicated,  the  single  ratios  yielding  the  distributions  of  the 
table  are  not  stochastically  independent  and  are, indeed,  highly  correlated, 
are  based  on  the  same  experimental  units  (convective  bands),  a»  sometimes 
result  from  stations  In  close  proximity.  Further  Investigation  is  needed. 
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either  in  analysis  of  Table  2,  or  throu!»h  confirmation  of  the  annarent  effects 
in  Table  2 by  more  ap^ronriate  statistical  methodoloRv. 

The  bi'C  reports  contain  interestinp,  'area  of  effect*  maps.  Tv;o  of  these 
maps  - one  for  composite  single  ratios  and  the  other  for  composite  double 
ratios  of  precipitations  are  shown  for  illustration  purposes  in  the  following 
two  pages,  "hen  a band  is  seeded,  silver  iodide  particles  'irlft,  diffuse, 
and  produce  ice  particles  which  in  turn  grow,  drift,  and  descend  to  earth. 

The  envelope  of  ground  interception  points  of  the  collection  of  particles  is 
termed  the  'area  of  effect  . The  heavy  dotted  line  forms  the  boundary  of 
the  expected  area  of  effect.  The  shaded  area  is  t'lat  subset  of  the  area  of 
effect  ^'herein  the  precipitation  increase  was  5n'<  or  more.  According  to  the 
reports,  extensive  portions  of  the  area  of  effect  bad  ratios  (both  S.’l.  & 
n.P.)  of  1.5  or  more  as  indicated  by  these  maos. 
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Figure  1:  Cooposite  doable  ratios  of  precipitation  for  1967-71,  all  bands, 

ground  seeded,  56  seeded  and  51  not-seeded. 


Figure  2:  Cooposite  double  ratios  of  precipitation  for  1967-71,  ail  bandt. 


This  section  would  be  incomplete  without  a summary  of  conclusions  arrived 
at  in  the  NWC  reports  [4,  7].  A brief  synopsis  of  the  main  conclusions  is  pre- 
sented below. 

Evidence  is  conclusive,  say  the  reports,  that  cloud  seedlnp  is  an  effec- 
tive means  of  increasing  precipitation.  Average  precipitations  due  to  seeding 
increased  by  50%  or  more  in  majority  of  target  area  stations.  There  was  an 
attempt  to  post-stratlfy  the  bands  by  their  alr-nass  sturctures  and  500  mb 
temperatures.  The  reports  conclude  that  average  precipitations  in  warm  and 
unstable  bands  increased  by  more  than  100%  and  covered  a much  larger  portion  of 
the  target  area.  Seeding  is  most  effective  when  500  mb  temperatures  are 
between  -17°C  and  -20°C.  There  is  some  evidence  to  the  effect  that  day-time 
seeding  is  less  effective  than  night-time  seeding.  *'ost  Importantly,  the 
reports  conclude  that  the  effect  of  seeding  is  through  an  Increase  of  hand 
duration  rather  than  an  Increase  of  band  intensity,  that  Is,  the  larger 
average  precipitations  in  seeded  bands  can  be  attributed  to  longer  band 
durations.  Data  summarization  in  this  report  and  in  ONR  Technical  Report 
No.  118  is  preliminary  to  reanalysls  of  the  Santa  Barbara  Data.  The  objectives 
of  the  reanalysis  are  both  to  evaluate  the  use  of  more  sophisticated  statis- 
tical methods  and  to  sharpen,  confirm  or  refute  conclusions  drawn  initially 
by  NAWC. 
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III.  DATA  PO”  RESPO’''SF  ‘IIWFACE  A”PI.TCATTf'*i'l 


We  review  data  used  below  In  this  report  in  annlicatlon  of  response 
surface  tnethodoloF.y . I'orth  American  ’leather  Consultants  had  tvro  projects 
(concurrently),  one  for  the  TJaval  Weapons  Center  at  China  Lake,  California, 
and  the  other  for  the  Pureau  of  Reclamation  at  Denver,  Colorado.  The 
analysis  undertaken  for  the  bureau  of  Reclamation  used  a >"uc!i  larger  network 
of  ralngage  stations  than  used  for  the  ilaval  'Weapons  Center  stud^'.  This 

Increase  in  the  number  of  ralneage  stations  resulted  from  an  extension  of 

the  target  area  to  the  east.  ’Thile  the  two  reports  discussed  in  Section  II 
were  done  on  behalf  of  the  Waval  ’’eapons  Center,  there  is  a third  report 
by  Bro',m,  Thompson,  and  Elliott  fl’,  which  was  done  for  the  Bureau  of  Re- 
clamation with  the  larger  network  of  ralngage  stations. 

The  first  of  the  tv'o  Waval  "eaoons  Center  reports,  written  in  February 
1D72,  had  only  80  raingage  stations.  A list  of  these  stations  is  contained 
in  Appendix  5 of  '’AT,  \11  single  and  double  ratios  in  this  report  were 

calculated  for  these  80  stations.  The  second  and  the  final  report  (October, 
1975)  had  112  ralngage  stations  listed  in  Anpendlx  A,  '7\  out  of  which  72 
were  common  with  the  list  in  A).  That  is,  from  the  list  of  stations  in  the 

former  report,  eight  stations  were  drooped  and  fortv  added  to  arrive  at  the 

list  in  the  latter  report.  However,  single  ratios  in  '^7''  are  reported  onlv 
for  100  stations,  that  is,  12  stations  listed  in  Appendix  A of  ’71  were 
again  dropped  from  the  analysis.  On  the  other  hand,  the  report  done  for  the 
Bureau  of  Reclamation  attempted  to  use  a network  with  about  178  ralngage 
stations,  100  of  which  were  common  with  '^P.'C  renort  ' 7 ' and  78  were  further 
additions.  The  reasons  for  dropping  or  adding  stations  are  not  clear,  '’^ere 
are  stations  which  have  been  dropped  from  the  analysis  but  have  precipitation 
data  available  for  them*  H9  and  MIO  are  examples  of  such  stations. 
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A sll<<htlv  different  Vlnd  of  nroblen,  potentially  capable  of  introducing  an 
element  (,f  nonco  p.. lability  vori'nis  :*nalv3es,  relates  to  r.lssln;'  obser-iratlons . 

There  are  a substantial  number  of  ralnp^ap.e  stations,  which  are  In  the  analyses, 
but  have  larpe  r t of  nlsslnp  nreclpltatlon  data.  In  fact,  there  are  few 
st.ntlons  which  have  precipitation  measurements  for  all  107  convective  bands. 
Clarifications  on  this  problem  were  souebt  from  "bV-’C  and  reasons  plven  were: 

(I)  some  ralnpapes  were  c>ianped  to  a different  tvoe  of  measuring,  de^'lce, 

(II)  some  ralnpapes  were  added  In  the  course  of  the  experiment,  and  (111) 
there  were  enulpment  failures. 

’“'V-'C  provided  a magnetic  tane  containing  data  on  all  variables  (excent 
location  coordinates  of  stations)  collected  during  the  course  of  t’-e  "'-’nta 
Barbara  experiments  and  used  either  for  M"C  analyses  or  for  the  Tureau  of 
Reclamation  analyses.  This  tape  consists  of  1"  files,  each  with  a different 
data  set.  Pile  has  precipitation  data  for  173  stations  from  ir>7  experi- 
mental units  (bands).  Locations  coordinates,  latitudes  and  longitudes  In 
degrees  and  minutes  and  altitudes  In  meters,  were  provided  by  for  each 

of  the  173  stations  separately. 

The  precipitation  and  location  data  first  considered  In  application  of 
the  response  surface  methodology  consisted  of  data  for  the  173  stations 
discussed  above.  The  locations  of  these  stations  are  summarized  In  Table  3. 


Table  3:  LOCATIONS  OF  RiUNOAOE  STATIONS  BY  LATITUDES  AND  LONGITUDES 


*Three  stations,  A45,  F6006,  and  E8832,  were  omitted  from  our  analyses  because  of  In- 
correct altitude  specifications. 


-16- 


Stations  labelled  Target  and  Control  In  Table  3 were  used  In  our  final  analyses 
with  these  two  areas  defined  by  the  Indicated  hounds  on  latitudes  and  lonf»ltudes* 
the  renalnlntr  stations  and  regions  were  omitted  because  it  was  lud-^ed  that 
orecipltatlon  measurements  were  too  sparse  or  that  they  "ere  too  far  from 
anticipated  regions  of  effect.  The  division  between  Target  and  Control  areas  re- 
m-^ins  a north-south  Up...  sTlehtly  to  the  vusf  I'f  the,  '!rourd-s««dln  rltv,  at 
latitude  34°32'  and  longitude  120*^01’  . Our  major  study  Is  based  on  107  Target 
area  stations  and  34  Control  area  stations. 

Appendix  Figure A-1  contains  a niap  showlne  locations  of  rany  of  the  ralngage 
stations  originally  considered  to  be  used  for  testing  the  response  surface 
approach  to  band  summarization  data.  Appendix  Tables  .A-l  and  A-2  give  a listing  of 
Tnfget  and  Control  area  stations  actuallv  used  together  with  their  three  location 
coordinates.  A similar  listing  of  stations  not  used  Is  contained  in  Table  A-3. 
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IV.  THE  CHOICE  OF  HESPOHRE  RHPFACER 


A.l  Introduction . Our  approach  to  the  sunmarl zatior  of  precinitation  data 
Is  in  tv'o  stages.  In  the  first  stape.,  a ’■esponsc  surface  is  fitted  to  the 
raingape  Tneasurements  for  each  convective  hand  separately  for  defined  control 
.and  target  regions.  The  Independent  variables  for  both  regions  are  simply 
the  raingage  location  coordinates,  latitude  and  longitude.  In  the  second 
stage,  integrals  yielding  volumes  under  the  surfaces  respresentative  of  total 
precipitations  for  specified  areas  of  interest  are  calculated  as  a sumnarlzatlon 
of  the  precipitation  data.  Thus,  various  calculated  precipitation  volumes 
may  be  derived  to  measure  the  response  to  a band  in  the  target  area  and 
similar  volumes  may  be  used  for  the  control  areas  as  possible  concomitant 
variables  describing  t>ie  strength’  of  the  band.  The  use  of  this  data 
sumnarlzatlon  will  be  reported  later  in  conjunction  with  concomitant  data 
summarized  bv  Hletson  in  OH"  Technical  Report  118. 

Preliminary  analyses  of  various  kinds  were  conducted  and  they  are  re- 
ported in  Sections  4.2  and  4.3.  These  preliminary  studies  were  based  on 
selected  bands  from  the  data  set  for  this  purpose.  Aspects  of  investigations 
of  goodness  of  fit  of  response  surfaces  "r  reported  also. 

Wot  all  possible  analyses  can  be  considered  and  inevitably  other  approaches 
would  be  selected  by  other  investigators.  Surface-fitting  techniques  are 
available  also  through  various  computer  packages  that  have  the  effect  of 
smoothing  responses  over  areas  of  Interest.  Such  a "local’  surface  fitting 
program  was  available  and  was  tried,  but  little  is  known  about  the  proba- 
bilistic Interpretations  possible  and  we  believe  that  these  routines  tend  to 
over-fit  data.  A discussion  of  such  methodology  as  used  In  m.eteorologv  Is 
given  by  Stephens  T ft  1 . The  use  of  these  procedures  was  rejected  but  thev  could 
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yield  summarization  methods  of  use.  The  basis  for  the  decision  was  that  high 
correlation  with  simple  averages  of  precipitation  were  anticipated  and.  Indeed, 
this  occurred  with  our  method  also. 

4.2  Pilot  Study  on  Selected  Banda.  The  appropriate  choice  of  regression 
equation  to  be  fitted  with  latitude  and  longitude  of  ralngage  stations  as 
Independent  variables  Is  Important.  A pilot  study  was  conducted  on  selected 
bands  to  determine  the  simplest  regression  equation  that  would  adequately 
approximate  data  from  all  bands. 

Bands  for  the  pilot  study  were  selected  to  be  representative  of  the 
totality  of  bands.  Twelve  bands  were  selected,  six  early  In  the  experimentation 
and  six  late.  Three  of  the  early  and  three  of  the  later  bands  were  seeded. 

The  convective  bands  used  for  the  pilot  study  were  bands  1,  2,  3,  4,  5,  7, 

94,  95,  96,  97,  98,  and  99.  The  six  seeded  bands  were  4,  5,  7,  94,  95,  and 
96.  In  the  pilot  study.  Station  N21  In  the  Control  Area  was  omitted  because  a 
corrected  latitude  specification  was  needed  and  later  supplied.  This  accounts 
for  minor  discrepancies  between  control-area  table  entries  for  the  pilot  study 
given  In  Tables  5,  A-4 , A-5  and  A-6  and  complete  summary  data  given  In  Tables 
A-8,  A-10  and  A-12. 

Early  work  showed  that  altitudes  of  ralngage  stations  were  not  helpful 
In  the  description  of  response  surfaces  nor  were  transformations  of  the 


precipitation  measurements.  Polynomial  regression  models  of  degrees  one, 
two,  and  three  respectively  in  latitude  and  longitude  were  fitted  to  ralngage 
measurements  for  all  the  twelve  bands  and  for  target  and  control  areas 


separately.  For  a given  band  and  region,  let  y^  be  precipitation  at  station 
J In  Inches,  be  latitude  of  station  j In  degrees,  Xj,^  be  longitude  of 

station  j In  degrees.  The  three  models  used  In  representation  of  y^  were 
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(3) 


"1  = ^no  ■"  ^lO^^ij-  ^ Yni(’'2i  “ ^^2^ 


'''  ^02^^21  *2^  ^ ’ 


(4) 


and 


^1  ' “on  “lo  ^’‘n  ■ “ni^’^21  ' *2^ 


- s2 


+ «2n^^]1  - ^ ■*■  “02^^2 


21 


- >2 


+ «12^^1  ■ ■ ^2^^  “03^’^2,i  ' ^2^^  ■"  "r 


(5) 


where  x^  and  x^  are  the  arltVimetic  means  of  latitudes  and  lonRltudes  of 
stations  In  use  for  any  particular  convective  hand  and.  In  each  model, 
is  a random  error.  For  each  model,  the  regression  coefficients  plven  by 
Greek  letters  must  be  estimated.  The  subscript  1 In  the  equations  above 
runs  over  the  n^  stations,  somewhat  different  from  band  to  band,  Xj^  and 
x^  are  also  different  for  different  bands  for  the  same  reason. 

The  pilot  studv  on  choice  of  model  was  Initiated  to  find  the  slnnlest 
adequate  model.  The  possibility  of  use  of  a model  of  bipher  deqree  than 
those  In  (3)  - (5)  was  open  and  It  had  been  honed  that  a second-deRree 
model  mlf'ht  be  adequate.  In  the  end,  model  (S)  was  necessary  and  use  of 
this  model  did  seem  suitable  in  that  data  summarization  was  our  Roal, 
the  cubic  model  fitted  substantially  better  than  the  quadratic  model, 
and  residual  variation  about  the  cubic  model  seemed  relatively  free  of 
systematic  effects. 
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Table  h f^lves  the  proeort  ton  of  precipitation  variation  explalne<' 
throuith  use  of  linear,  nuadratic,  and  cubic  response  surfaces  in  the  target 
area  for  each  of  the  tvelve  test  bands.  This  table  also  renorts  on  the 
mean  precipitation  v and  the  measure  of  unexplai.ied  variation  after  fitting 
cubic  response  surfaces  (”S'^'-cubi c) . Similar  information  for  regressions  in 
the  control  area  is  "iven  in  Table  5.  Computations  of  the  various  nuantltles 
reported  in  Tables  and  5 have  been  done  using  standard  regression  methods. 
Equations  (3)  - (5)  verc  fitted  as  shovm  because  initial  failure  to  center 
the  Independent  variables  and  x^  led  to  computational  problems.  Pradlev 
and  Srlvastava  reported  on  this  matter  in  nCT  Technical  "eport  111. 


Table  A:  SlIfP'APY  STATTr>TTC'l  fOH  PnLY’'’ovTAL  Rrr.HEfiSTn’ir,  ON 

T ■ TAPGET  area  p'CR  EELECTEt)  PA'IOS 


Band  .’o.  n 

y 

2 

Lr^EAR 

nUADRATTC 

' i 

CIIRTC 

’ffio  (niRic) 

1 1 

75 

.0637 

.00484 

.21377 

. 34970 

.45271 

.00302 

2 

75 

.1029 

.00781 

.2034? 

.37915 

.49288 

j .00451 

3 

75 

. 1665 

.00709 

.13967 

.15379 

.00949 

t* 

75 

.0376 

.00116 

.04631 

1 

.20771 

.26880 

.00096 

5* 

RA 

.15’’7 

.00323 

.16722 

.2690? 

.56887 

! .00397 

7* 

8C 

.035? 

.00259 

.32275 

.39495 

. 50381 

.00142 

94* 

91 

.0170 

,nno^2 

.31601 

.48785 

.58231 

.00029 

95"  ' 

')/, 

. ] 320 

.0129  ■ 

.45663 

.52864 

.59322 

' .00580 

96* 

03 

.6634 

.17825 

.58399 

I 

.65960 

.75583 

.04025 

97 

93 

.4362 

.13965 

.4831  2 

.59804 

.64237 

.05536 

9f. 

93 

.0940 

00769 

.20961 

.31405 

.36529 

.00542 

qo 

92 

.•325E 

I 

C 

c 

.28664 

.37727 

.45371 

t .10916 

*Hands  with  asterlsbs  are  seeded. 
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Table  5-  SITCIAOV  STATISTICS  FOR  POCYirOMIAL  RECRFSSIO'’S  ON 


’ CO  'trnL  ArCA ' 

'O^  SELECTEO 

BA.NnS 

Band  No. 

n 

V 

s2 

LINEAR 

OUADRATIC 

COBIC 

MSE  (CUBIC) 

1 

15 

.1013 

1 

010^0  ^ 

.31404 

.67735 

.74744 

.00771 

2 

15 

.2540 

00776 

.42910 

. 51983 

.70191 

.00647 

3 

16 

. 0900 

00331 

.22049 

.63129 

.91370 

.00071 

4* 

16 

.0563 

0012O 

.01738 

.46845 

.88551 

. 00034 

5* 

16 

.1450 

00337 

1 

.53032 

.80450 

.93254 

.00141 

1* 

16 

. 1125 

004  30 

.43388 

.79435 

.36047 

.0016S 

94* 

28 

.0693 

00342 

.24475 

.29054 

.35515 

.00331 

95* 

32 

.1138 

01290 

. 37361 

. 60504 

. 60594 

.00717 

9b* 

32 

. 3009 

05700 

.62736 

.74825 

.79654 

.01634 

97 

32 

.2822 

04215 

. 36790 

.55173 

.58024 

.02493 

98 

32 

.0622 

00162 

.05370 

. 16033 

.23034 

.00164 

99 

31 

.6835 

03151 

.03765 

. 17507 

.26793 

.03294 

*Bands  with  asterisks  are  seeded. 


Values  of  P in  Tables  and  5 measure  the  fraction  of  variation  in 

precipitation  measurcnents  explained  bv  linear,  ouadratlc  and 

cubic  models  for  each  band  for  target  and  control  areas.  In  general, 

substantial  improvement  resulted  for  use  of  the  ouadratlc  model  over  the 

linear  and  for  use  of  the  cubic  over  the  quadratic.  The  decision  to 

use  the  cubic  mode;]  rather  than  to  proceed  t«'  a higher  degree  model  came 

also  from  examination  of  residual  variation  as  discussed  below.  Indications 

of  heterogeneous  error  variability  in  both  the  observed  data  and  the 

2 

residual  errors  comes  respectively  from  the  values  of  s and  mre  (Cubic) 
columns.  This  indicated  heterogeneous  varlabilltv  must  he  considered 
in  further  analyses  and  suggests  use  of  wel'd’ts,  for  example,  in  planned 
covariance  analyses. 
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riottlnf,  rocliniqut's  v'cre  used  extenslvelv  In  determination  of  choice 
of  response  surfaces.  It  is  not  possible  to  include  comnlete  sets  of 
plots  hut  results  are  illustrated  for  Rand  Tarpet  area,  a seeded 
band  with  relatively  hlph  precipitation,  Appendix  Fipure  A»2  '’ivos  a plot 
of  the  precipitation  tgta . Unear,  quadratic,  and  cubic  response  surfaces 
as  fitted  to  the  data  ereriven  In  Anpendix  Flp,ures  ''v-3,  an''  ‘.-5.  In  Figure 

VS,  it  is  apparent  that  the  cubic  surface  becomes  negative  in  the  lov/er  left 
comer  of  t^  tareet  ropion.’-'ut  recall  that  this  portion  of  the  region 
is  not  a land  are.i  and  of  no  consequence  in  further  use.  All  these 
surfaces  are  oi  i'  r over  regions  of  laior  interest.  Figures  A-6,  A-7,  and 
A-v,exhlMt  r-sldual  var  lat  ion,  lepari  ures  of  observed  precipitations 
from  the  fitted  surfaces.  In  general  over  the  test  bands,  linear  and 
omdratii  models  tended  to  leave  residuals  with  some  svstc^natlc  patterns, 
both  In  re  .ird  il.'.ns  ire'  -"'nltudes.  Residuals  from  cubic  surfaces 

'/•■re  rel.aflvel.  at  net  <onpletelv  free  of  such  systematic  effects. 

Som«‘  furt'-.  t . yi'  !-  itl  >n  : r>sldiials  fol  lo^.'s  In  the  next  subsection. 

4.J  Pr<- 1 1 V 1 n,,r  .•  Istrihutions  and  Transformations. 

This  subs.  . f , ! . . : • ii  .-x.imlnat  ion  of  distributional  properties  of 

(1)  prerlpifati  - i-  ireiaer.f  , n r>^-;s  stations  within  a band  and  (il) 
resicual  l.it  1 tr  m liu-ar  qua.lritli  tnd  i.ubic  response  surfaces.  In 
addition,  we  teg'  rt  'rleflv  on  the  use  of  transformations. 

Two  types  ol  studies  th.it  directly  consider  distributions  of  raingage 
measurements  exist.  f>n  type  of  study  fits  specific  distributions, 
usu.illy  CansTLi  or  log  normal,  to  precipitation  data.  A consequence  of 
such  study  'nay  he  to  select  a supposedly  appropriate  statistical  test, 
perhaps  a nonparanetrlc  one  to  examine  seeding,  effects.  Examples  of  this 
can  be  found  in  hennis  and  Schock  f2'  and  in  Ouran  and  Mlelke  TBl.  Hanson 
in  OMF  Technical  Report  110  has  considered  asymptotic  comparisons  of 
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efficiencies  of  certaii\  rank  tests  for  this  project.  \ second  tvpe  of 
study  considers  transformations  to  "normalize  data.  '^mlth,  \dderly  and 
Botln.^alte  report  on  use  of  the  square-root  transformation  In  an 

exporlnent  In  South  Australia. 

We  applied  two  transformations  to  precipitation  data  to  determine 
if  these  data  .showed  anv  evidence  of  beconlnp.  more  nearly  normally 
distributed  followln.?  transformation  than  before.  The  first  transfonr.atlon 
Is  loo.arl thnl c and  commonly  used.  The  second  is  a double  lor»arlthnic 
transformation  devised  after  the  first  was  judpcd  unsatisfactory.  The  exact 
‘'or-  s are 


X = In  (1  + Y)  (6) 

and 

Z = In  (1  + X),  (7) 

where  Y Is  the  observed  precipitation.  Tables  A-A  to  A-b  ''hor;.  for  tl  e fclve 
test  bands  and  control  and  tarpet  areas»  values  of  the  variance, 
standardized  skcvn,>ss .and  standardized.  Juirto-il'  for  the  precipitation  data, 
the  z - transformations , see  (7),  of  the  precipitation  data,  and  the  residuals 
from  the  cubic  surface  fitted  to  precipitation  data.  I''ata  is  on  file 
in  similar  ce  rations  for  the  transform  (b)  and  for  residuals  from  linear 
and  ciuadratlc  surfaces  results  are  not  reported  hero  because  they  are 
Intermediate  to  those  shown. 

The  asvmmctrlc  and  usually  hlph  peakedness  of  precipitation  distri- 
butions motivated  tlie  choice  of  transformations  (b)  and  (7).  It  Is  knovm 
that  loKarithmlc  tran.sform-it Ions  Iielp  to  normalize  skev'ed  and  leptoUurtlc 
distributions.  Hov/ever,  In  our  case,  transformatlon.s  chosen  helped 
only  marginally  to  shift  downwards  the  hlglier  moments  of  the  distributions. 
The  drops  In  variance  and  skewness  are  too  small  to  be  of  any  consequence 
and  tile  reductions  in  kurtosls. even  though  substantial  for  some  bands,  are 
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far  from  enough  Lo  make  then  close  to  those  for  a normal  distribution. 

The  conclusion  is  that  the  transformed  variables  and  7.  are  not  normal. 

These  results  arc  consistent  with  those  reported  in  the  llter-’turo.  '■  further 
conclusion  is  that  any  benefit  resulting,  from  transformation  before 
surface  flttinp  is  more  than  offset  by  the  interoretdtive  difficulty 
that  would  result  if  surfaces  were  fitted  to  transformed  precipitations. 

Some  other  transformations  vrere  tried  also  but  were  not  helpful. 

Investigations  into  the  distributional  aspects  of  residuals  for 
stations  within  a band  from  fitted  rt^sponse  surfaces  were  also  undertaken. 

The  residual  at  station  j is  ilefined  as  - y^  where  y^  is  the 
observed  precipitation  at  station  j and  v^  its  estimated  value  using 
the  appropriate  response  surface.  It  is  clear  from  Table  A-G  that  the 
residual  distributions  in  the  target  area  are  always  right-skewed.  The 
kurtosis  for  residuals  from  all  bands  and  all  response  surfaces  were 
also  positive,  and  sometimes  very  high  in  magnitude.  One  conclusion  is 
that  residuals  from  cubic  fits  had  right-skewed  and  extremely 
leptokurtic  distributions  in  the  target  area. 

Ue  have  made  the  decision  to  fit  nreclpltatlon  responses  without 
transformation  to  response  surfaces  which  are  the  general  cubic  in 
two  dimensions,  latitude  and  longitude.  Other  decisions  could  have 
been  made  and  other  investigations  could  have  been  conducted.  At  some 
point,  decisions  must  be  made.  Justl^^lcatlon  of  our  decision  is  that 
our  primary  purpose  in  the  fitting  of  response  surfaces  was  for  data 
summarization.  The  option  of  transforming  response  statistics  calculated 
from  the  surfaces  is  still  open  and  will  be  considered. 

4.4  Correlations  in  Polynomial  Regression.  It  was  detected  early  in 
the  surface  fitting  v;ork  that  the  Independent  variables  (latitude  and 
longitude)  vfere  such  that  the  absolute  sample  correlation  coefficients 
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between  each  of  them  and  their  powers  turned  out  to  be  close  to  unity. 

This  caused  difficulty  in  obtalninR  least  squares  estimates  of  regression 
coefficients.  The  problem  was  investigated  by  Hradley  and  Srivastava  and 
they  have  reported  in  OitP  Technical  Report  111.  '.^1110  the  problem  has 

been  noted  by  others,  it"  rationalization  may  be  of  value  both  in  research 
and  teactiing.  The  difficulty  v;as  removed  on  this  project  through  centering 
of  the  Independent  variables  in  the  computational  'rorV. 
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V.  CUI-'IC  RESPONSE  SURE^CES  A 'ri  ESTIMATION  OF 
VOLIP'ES  OF  PRECIPITATION 

5.1  Introduction . Since  the  pilot  study  sMo\’ed  that  cubic  response 
surf.ices  provide  adequate  description  of  ralnpape  neasurenerits  for  a 
given  region  and  band,  surface  fitting  work  was  completed  for  all 
experimental  units  (bands)  and  target  and  control  regions  '’1th  use  of 
two-dimensional  cubic  response  surfaces  with  latitude  and  longitude  as 
Independent  variables.  These  fitted  surfaces  <jere  then  integrated  over 
designated  geographical  areas  to  evaluate  the  volume  of  rainfall  in  that 
area  for  tliat  band.  This  section  first  summarizes  the  cubic  surface 
fitting  'work,  and  then  goes  on  to  discuss  evaluation  procedures  for 
rainfall  volumes  together  vrlth  their  variances. 

5.2  Cubic  Surfaces . The  cubic  response  surfaces  with  latitude  and 
longitude  as  independent  variables,  were  fitted  to  all  107  experimental 
units  (bands)  and  for  target  and  control  areas  separately.  The  poly- 
nomial regression  model  used  to  fit  cubic  surfaces  is  given  by  the  equation. 
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for  observation  (station)  1,  where  y^ , , and  X2^  are  precipitation 

measurement,  latitude,  and  longitude  at  station  J,  j “ 1,  2,  ...,  n. 


and  Cj  is  the  random  error  associated  with  the  response  y^ . The  constants, 
Oj  and  ^2,  represent  the  centetlng  locations  used  for  latitude  and 
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lonKltiidc.  The  values  chosen  for  and  vrorc  34.0°  and  119.0°  for 
t.argot  area  surfaces  and  34.7  and  120.3  for  control  area  surfaces. 

The  estimates  of  rer.ression  coefficients  6^^  (0<£<3,  0<k<3-£) 
v;ere  obtained  usln>;  least  square  procedures. 

Note  that  (8)  differs  from  (5)  In  that  fixed  centering  at  and  near 
the  moans  of  the  corresuondinf;  variables,  Is  used.  This  Is  more  convenient 
for  tabluatlon  and  future  use. 

The  two  main  statistics  obtained  to  evaluate  the  Roodness  of  fit  of 

cubic  response  surfaces  ’-ere  t'^e  rercentaRe  of  variation  In  original 

measurements  explained  by  the  fitted  surface  and  the  residual  mean 

square  error  (USE),  the  estimate  of  variability  about  the  fitted  surface. 

Table  A-7  gives  for  target  area  surfaces  values  of  (1)  the  number  of 

stations  In  the  band , (H)  the  precipitation  average  for  stations  operative 

2 

In  the  target  region  for  the  band,  (ill)  the  sample  variance,  (iv)  100  R , 
the  percentage  variation  explained  by  the  fitted  surface,  (v)  the  MSE 

•k 

and  (vi)  the  calculated  value  of  the  F-ratio  . The  corresponding  infor- 
mation for  fitted  surfaces  In  the  control  area  is  contained  In  Table 
A-8.  The  estimated  regression  coefficients  b^^^  are  listed  bv  bands  in 

Table  A-9  for  the  target  area  and  In  Table  A-10  for  the  control  area. 

2 

Consider  values  of  lOOR  for  the  target  area  In  Tahlo  A-’?-. 


Under  classical  assumptions  for  regression  models,  F has  the  Snedecor 
variance-ratio  distribution.  V.lille  significance  levels  are  referenced 
In  the.sc  tables,  validity  of  use  Is  In  question. 


They  broadly  confirm  the  pattern  observed  in  preliminary  study  of  the 
test  bands.  The  percentage  of  precipitation  variation  explained  varies 
from  a low  of  13.8%  for  band  67  to  a high  of  84.1%  attained  by  bands 
11  and  12.  Tliere  are  as  many  as  32  bands  for  which  the  fitted  surfaces 
did  not  account  for  even  40%  of  the  variation  in  the  original  observations. 
On  the  other  hand,  30  hands  had  values  of  100^2  excess  of  60%.  The 
performance  of  cubic  response  surfaces,  as  judged  by  the  values  of  lOOR  ’ , 
can  be  considered  to  be  fairly  good. 

2 

It  was  conjectured  that  values  of  R and  MSE  may  be  dependent  on 

the  magnitude  of  average  hand  precipitation.  There  was  no  relationship 

2 

whatsoever  between  R and  the  simple  average  y,  but  a strong  relationship 
between  y and  !!SE  was  evident. 

The  values  of  100" 2 r^r  cubic  surfaces  in  the  control  area  are  nearly 
all  very  high.  Thus,  data  description  by  means  of  cubic  response 
surfaces  in  the  control  area  fitted  well.  This  is  reflected  also  in 
unusually  small  and  almost  constant  values  of  MSK's  (when  rounded  to 
two  decimals  places).  But  there  are  fewer  raingage  stations  in  the 
control  area  and  some  of  the  apparent  goodness  of  fit  is  due  to  the 
large  number  of  parameters  estimated  relative  to  the  number  of  obser- 
vations available.  Again,  there  was  a strong  dependence  of  MSE's  on 
the  y's. 

5.3  Volumes  of  Precipitation  in  Designated  Geographical  Areas.  The 
objective  of  the  response  surface  work  was  the  sunmarlzation  of  band 
precipitation  data.  The  summary  measures  to  be  used  are  the  volumes  of 
precipitation  under  the  fitted  response  surfaces  for  designated 
geographical  areas  of  Interest. 
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We  Illustrate  precipitation  volume  calculations  through  use  of  (8)  and  for 
rectangular  areas  with  north-south,  east-west  orientations.  Other  areas  may  be 
used  with  only  minor  complications  and  may  be  used  later.  It  is  assumed  that  any 
area  used  is  completely  in  either  the  target  or  control  region  and  that  the  appro- 
priate estimated  surface  from  (8)  Is  used.  The  estimated  surface  may  be  written. 


y(Xi,  x^)  = (x^-a^)  + b^^  (x2-a2>  + ^>20 


2 3 

bfi  (x^-a^)  (x^-a^)  + b^2 


(9) 


2 2 3 

^21  (x2-a2>  + bj^2  (x2-a2>  + b^^^  (x^-a^)  • 


Let  us  assume  that  the  precipitation  volume  V in  Inches  of  precipitation 
per  degree  of  latitude  by  degree  of  longitude  Is  required  for  the  region, 
a < Xj  < b,  c < X2  < d or  for  a < B.  y < *2~“2 

ft  «»  b-Oj^,  y = c-a2,  9 “ d-a2.  The  precipitation  volume  is  obtained  by  elemen- 

tary integration  as 

V = bg^  (B-a)(S~y)  + b^g  (ft^-a^) (6-y) /2  + bg^  (ft-a) (6^-y^) /2  (10) 

+ b2g  (6^-a^)(6-Y)/3  + b^^  (ft^-a^) (6^-Y^) /4  + bg2  (6-a) (6^-Y^ /3 

+ bgg  (ft^-a^)(6-Y)/4  + b2j  (B^-a^)(6^-Y^)/6  + b^2  (6^-Y^ /6 

+ bgj  (6-a) (6^-Y^)/4, 

a linear  function  of  the  estimated  regression  coefficients.  But  the  variances  and 
covariances  of  these  regression  coefficients  may  be  estimated  and  hence  the  esti- 
mated variance  of  V may  be  obtained. 

Volumes  and  their  variances  were  evaluated  for  each  band  for  the  following 
areas'  (1)  the  complete  tcrget  area,  34.0  < Xj^  < 35.25,  118.0  < X2  < 120.02, 

(11)  that  portion  of  the  target  area  in  Santa  Barbara  County,  34.4  < Xj  < 35.0, 
119.51  < *2  ^ 120.02,  (ill)  that  portion  of  the  target  area  which  enclosed 
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stations  in  most  of  Ventura  and  most  of  Los  Angeles  counties  with 
34-0  < < 35.0,  118.0  ^2  119.51,  (Iv)  the  sum  of  areas  defined 

In  (11)  and  (111),  and  (v)  the  control  area,  34.4  < < 35.25, 

120.02  ^2  ^ 120.60.  The  estimated  volumes  and  their  variances  for 

the  above  five  areas  are  presented  In  Table  A-11.  It  should 
be  pointed  out  here  that  tarf»et-area  response  surfaces  were  used  for 
evaluations  of  the  first  four  volumes,  control-area  surfaces  were  used 
for  evaluations  of  volumes  (v) . 

Estimates  of  volumes  bv  bands  were  compared  with  simple  averages^of 
ralngage  measurements  from  stations  In  the  designated  areas  of  Interest. 
Table  6 shows  values  of  the  correlation 


Table  6:  CORRELATION  COEFFICIENTS  BETWEEN  VOLUMES 

AND  SIMPLE  AVERAGES  OF  RA INGAGE  MEASUREMENTS. 


each  correlation  coefficient  in  the  table  should  be  baaed  on  107  volume- 
average  pairs  corresponding  to  107  hands.  However,  one  band  In  the 


These  averages  and  their  variances  are  given  In  Table  A-12. 
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target  area  and  six  bands  In  the  control  area  had  negative  volumes  and 
hence  were  deleted  from  the  calculations.  The  negative  parts  of  surfaces 
used  were  off  the  land  area  and  redefinition  of  these  regions  will  be 
necessary . 

Table  6 shows  that  all  correlation  coefficients  are  very  high.  The 
implication  of  these  high  correlations  is  that  tlie  simple  average  of 
raingage  measurements  may  be  an  adequate  summary  measure.  Further 
analysis  is  likely  to  involve  v/elghtlng  of  summary  measures  through  use 
of  their  variances,  whether  or  not  there  is  advantage  in  use  of  precipi- 
tation volumes  over  simple  averages  may  depend  on  the  effects  of  weighting. 
Both  types  of  measures  will  be  considered  in  examination  of  the  effects 
of  cloud  seeding.  The  investigators  did  judge  that  it  was  important  to 
consider  the  possibility  of  more  informative  data  summarlzatton  than  Che 
use  of  simple  averages. 
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VI  SUMMARY  AND  CONCLUSIONS 


This  report  is  on  one  of  a series  of  investigations  on  surnmarization 
of  data  in  weather  modification  experiments  with  particular  reference  to 
the  Santa  Barbara  experiment.  Summary  measures  of  precipitation  have  Peen 
constructed  using  a two-stage  procedure.  In  the  first  stage,  response 
surfaces  v;ere  fitted  to  ralngage  measurements  for  each  experimental  unit 
with  ralngage  location  coordinates  as  Independent  variables.  In  the 
second  stage  Integrals  of  these  surfaces  were  evaluated  to  obtain  measures 
of  total  rainfall  over  designated  areas  for  each  experimental  unit. 

The  method  proposed  is  tested  with  data  made  available  by  Iiorth 
American  Weather  Consultants  from  Phase  I of  their  Santa  Barbara  Con- 
vective Test  Seeding  Program.  A review  of  these  experiments  and  of 
reports  [1,  4,  71  based  on  data  generated  from  these  experiments  Is 
contained  in  Sections  n ^nd  II J. 

A pilot  study  with  twelve  selected  bands  was  undertaken  to  determine 
the  advantages  and  disadvantages  of  different  functional  forms  of  response 
surfaces.  Details  of  this  pilot  study  are  reported  in  Section  IV.  It  has 
been  shoT-m  (Section  4.2)  that  appropriate  response  surfaces  for  Santa 
Barbara  data  are  third  degree  polynomials  in  latitude  and  longitude  of 
ralngage  locations.  The  cubic  response  surfaces  did  produce  mean  squared 
errors  (HSE's)  which  were  substantially  smaller  than  the  sample  variances 
of  the  original  observations  (See  Tables  4 and  5).  The  patterns  of 
residuals  from  fitted  surfaces  were  examined.  It  has  been  shown  in 
Section  4.3  that  distributions  of  residuals  were  non-normal.  These  dis- 
tributions were  usually  right-skewed  and  leptokurtlc.  Section  4.3  also 
contains  discussion  on  distributions  of  the  original  precipitation  observations 
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and  of  two  transforms  of  then.  The  precipitation  observations  had  rl(?,ht-skewed 
and  leptokurtlc  distributions.  The  two  transformations  used  p.ave  only  Tnarplnal 
improvement  In  the  shapes  of  distributions. 

The  cubic  response  surface  fitting  was  extended  to  all  the  107 
experimental  units  and  for  target  and  control  areas  separately.  These 
surfaces  were  Integrated  over  designated  areas  of  Interest  to  produce 
total  volumes  of  rainfall  for  those  areas.  The  geographical  areas  consi- 
dered were  total  target  area,  three  subsets  of  the  target  area,  and  all 
of  the  control  area.  The  variances  of  the  computed  volumes  for  each  of 
these  areas  were  obtained.  Discussion  on  surface  fitting  and  evaluation 
of  volumes  is  given  in  .Section  V-  tables  of  results  are  given  In  A-7 
to  A- 11. 

The  percentages  of  variation  explained  by  cubic  response  surfaces 
were  moderately  good  In  the  target  area  and  were  excellent  In  the  control 
area.  Estimates  of  residual  variances  produced  were  substantially  lower 
than  the  variances  In  the  original  data.  Finally,  estimates  of  volumes 
for  selected  geographical  areas  were  found  to  be  highly  correlated  with 
simple  averages  of  ralngage  measurements  In  those  areas.  Both  precipi- 
tation volumes  and  area  averages  will  be  used  In  further  Investigation. 

Primary  use  will  be  In  association  v;lth  concomitant  variables  in  the 
evaluation  of  the  effects  of  cloud  seeding. 
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FIGURE  A-1:  PRO.IECT  >UP  SHOWING  KAINGAGI,  lOCATIONS,  IGVOAR,  AND  SEEDING  SITE 

SOURCE:  r.iKf  (3). 


FIGURE  A-3:  GRAPH  OF  LINEAR  RESPONSE  SURFACE:  BAND  96,  TARGET  AREA. 

(VERTICAL  AXIS  MEASURES  2.3  TIMES  THE  PREDICTED  PRECIPITATION  IN  INCHES). 


FIGURE  A-4;  GRAPH  OF  QUADRATIC  RESPONSE  SURFACE:  BAND  96,  TARGET  AREA. 

(VERTICAL  AXIS  MEASURES  2.3  TIMES  THE  PREDICTED  PRECIPITATION  IN  INCHES). 
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FIGURE  A-5:  GRAPH  OF  CUBIC  RESPONSE  SURFACE:  BAND  96,  TARGET  AREA. 

(VERTICAL  AXIS  MEASURES  2.3  TIMES  THE  PREDICTED  PRECIPITATION  IN  INCHES). 
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FIGURE  A-7 

GEOGRAPHICAL  DISTRIBUTION  OF  RESIDUALS  FROM  QUADRATIC 
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FIGURE  A-8 

OEOGRArHICAL  DISTRIBUTION  OF  RESIDUALS  FROM  CUBIC 
RESPONSE  SURFACE:  BAND  96.  TARGET  AREA 

(THE  VERTICAL  SCALE  IS  IN  INCIiES  AND  SHOWS  DOUBLE  THE  ACTUAL  RESIDUAL). 
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TABLE 

A-3;  INDEX  TO  PREC 1 UTATION  SIATIONS  NOT  INCI-UDED 
IN  nin  STUDY 
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TABLE  A-4:  VARIANCE  AND  SH^PE  STATISTICS 

EROM  SELECTED  BANDS 


Target  Area 


Band  No. 

Variance 

Skewness 

kurtosi s 

1 

.005 

1.850 

4.458 

2 

.008 

1.912 

4.779 

3 

.010 

2.251 

7.123 

4 

.001 

1.362 

2.563 

5 

.008 

0.390 

0.768 

7 

.003 

2.048 

4.845 

94 

.001 

1.472 

1.290 

95 

.015 

0.715 

0.047 

96 

.173 

-0.031 

-0.935 

97 

.140 

0.687 

-0.657 

98 

.008 

1.54S 

2.678 

99 

. 180 

0.770 

1.163 

_.J 
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FOR  PRIXIPI'^ATTONS 


CoAtrol  Area 


Variance 

Skewness 

Kurtosis 

.011 

0.801 

-0.915 

.008 

C.043 

-0.992 

.003 

0.621 

-0.250 

.001 

-0. 171 

-0.685 

.008 

0. 108 

-0.720 

. 005 

0.818 

0.675 

.003 

0.420 

-0.854 

.013 

2.085 

5.099 

.057 

1.267 

1.996 

.042 

1.799 

3.989 

.002 

0.321 

-0.807 

.052 

1.483 

3.466 

TABLf-:  A-S;  VARIAN'G'  /vNl)  SHAPf  STATISl’ICS  FOR  1.  - logll  + log(l  + Y)] 
FOR  SFLCrini)  BANOi; 


Band  No. 

j larger  Area 

, V.irianLt’  Oktwncss 

Ku.'-to.sis 

Control  Area 
Variance  Skewne.ss 

Kurtosis 

I 

1 . 00^ 

1 . 387 

2.243 

j . 00/ 

0.690 

-1.0G3 

n 

A. 

.004 

1.362 

2.131  * 

1 . 003 

-C.147 

-0.938 

.■5 

.0'J4 

1.482 

3.433 

1 .002 

0.431 

-0.531 

4 

.001 

1.126 

1 . 658  1 

1 .001 

-0.287 

-0.716 

5 

.005 

0.093 

-0.053 

i .005 

-0.163 

-0.772 

7 

.002 

1.725 

3.  i47 

.003 

0.480 

0.232 

94 

.001 

1 . .396 

0.939 

.003 

0.277 

-1.042 

95 

. 007 

0..358 

-0.726 

.006 

1.459 

2.316 

96 

.038 

-0.664 

-0.729 

.018 

0.454 

-0.330 

97 

.035 

0.150 

- 1 . 1 '4 

.012 

0.864 

1.027 

98 

.005 

1.102 

1.137 

.001 

0.203 

-0.852 

99 

I .022 

I 

-0.481 

0.362 

.004 

0.794 

1.145 
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TABLE  A-6:  VARIANCE  AND  SHAPE  STATISTICS  FOR  RESIDUALS  FROM  CUBIC 


SURFACES  FOR  SELECTED  BANDS 


Band  No. 

Variance 

Target  Area 
Skewness 

Kurtosis 

Variance 

Control 

Skewness 

Area 

Kurtosis 

1 

.003 

1.680 

5.392 

.003 

-0.581 

0.783 

2 

.004 

1.286 

5.461 

.002 

-0,540 

0.561 

3 

.008 

2.312 

7.457 

.000 

0.216 

0.734 

4 

.001 

1.584 

4.559 

.000 

1.367 

2.187 

5 

.004 

0.630 

3.324 

.001 

0.137 

-0.904 

7 

.001 

1.619 

6.898 

.001 

-0.312 

-0.713 

94 

.000 

0.872 

1.355 

.002 

0.279 

0.539 

95 

.005 

1.478 

4.422 

.005 

3.089 

13.348 

96 

.044 

0.864 

1.889 

.012 

0.811 

0.136 

97 

.050 

0.496 

1,058 

.018 

2.601 

9.755 

98 

.005 

1.467 

2.419 

1 .001 

0.889 

0.313 

99 

.098 

0.798 

1.742 

1 

1 .023 

1 

1 

1.123 

1.795 
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TABU-;  A-7;  PRECIPITATION  /CND  CUBIC  RESPONSl  SURPACt  STATISTICS 


TAkCrr  ARP A 


Band  No. 

No.  of 
Stat ions ,n 

Me  an 

Prccip. y 

Precip. 

Var. 

lOOR^ 

USE 

1 

75 

05^7 

.0046 

4 5.  2 

.00332 

5.9  o54 

75 

lu  26 

. Ou  76 

3 

.GC45  1 

7.0246 

d 

75 

1655 

. OCSS 

1 5.  5 

.00949 

1.3210 

4 

75 

j376 

. Q C1.2 

? 7.  1 

.00096 

2.6787 

5 

64 

15?/ 

. U !l  a 2 

5 7.  0 

. 0 0 .3  J 7 

10.887/ 

6 

84 

. 0 39? 

.0027 

41.2 

.00181 

5.7678 

7 

68 

.03  5? 

.0026 

50.9 

.001 42 

8.9703 

0 

e 3 

. 00  2? 

. 0 00  1 

36.7 

.00004 

4.7014 

9 

83 

.3075 

. 0 00  3 

6 3.6 

. 000  14 

14.1623 

1 0 

8 3 

.3146 

. 3 Ou  3 

5 4.  .1 

.00017 

9.6  416 

1 1 

8 

. 11  34 

.0  130 

84.  1 

.00231 

42.9631 

12 

8 3 

.14  32 

• 0 lo  3 

8 4.  1 

.00291 

42.7867 

1 3 

83 

. 0 3 5 a 

.0018 

64.  0 

.00074 

14. 4431 

1 4 

8? 

. 0 2 98 

. 0 0 26 

68.  7 

.00091 

l'^  .5228 

15 

82 

. 36o2 

. 0'"44 

26.  6 

.0  4 3 71 

3.1974 

16 

8 2 

. 0857 

. u 0 fi  6 

34.  4 

. 00o52 

4,1959 

1 7 

65 

.’4  33 

. 0 56  6 

7 8.  9 

.01337 

31  .2493 

1 3 

74 

. 3606 

. CCu3 

5 5.7 

.on  o 16 

9.6238 

19 

8 1 

. 5642 

. 1 054 

35.  1 

.077  y 3 

4.2697 

20 

8 1 

. 0336 

. 0 023 

51.6 

.00123 

9.4164 

21 

81 

. 0868 

. 0 0 7 1 

14.2 

.006  87 

1.3060 

22 

81 

.0,36  7 

.0  016 

40.8 

.00110 

5 .4312 

2 3 

69 

. 0954 

. 0041 

72.  1 

.00132 

16.9629 

24 

C9 

.1094 

. 0 046 

6b.  1 

.001 76 

12.7843 

25 

. C547 

.6017 

5 7.  1 

.00085 

8.5667 

?6 

69 

. 05  77 

. J 0 3 1 

36.  7 

.00223 

.3.8  05  7 

27 

75 

. 0091 

. C C02 

27.  8 

.00019 

2. 7760 

28 

75 

. 1267 

. U 1 0 4 

50.6 

.f’05  86 

7,4044 

?o 

■'9 

. 0u94 

. 0 0 07 

6 5.  2 

.00026 

1 4 . 3 49  7 

30 

81 

.1625 

. 0 04  3 

5 6.5 

.00456 

10.2636 

Footnote:  Vsigni ficance  probabilities  of  F-ratios  can  be  approximated 

fror.1  the  standard  F-table  with  9 and  (n-10)  degrees  of  freedom 
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I TABLE  A-7:  PRECIPITATION  AND  CUBIC  RESPONSE  SURFACE  STATISTICS, 

TARGET  AREA  (CONTINUED) 


Band  No. 

No.  of 
Stations 

Mean 

,n  Prccip.y 

Precip. 
Var. s ^ 

lOOk^ 

MSE 

31 

.0127 

. 0005 

1 7.8 

. 000  4.3 

1.7341 

32 

84 

.0224 

. Cull 

35.  8 

• oon  77 

4.5846 

33 

88 

. 0522 

. 0 034 

b2.  5 

.00142 

14.4713 

3A 

d's 

.3091 

. 0235 

6 4.0 

.00928 

iP'.jieo 

35 

89 

. 23  6U 

.01-^2 

4 4.8 

.01059 

7.1154 

36 

89 

. 0t)4C 

.0  j40 

49.  7 

.00225 

8 .6699 

37 

92 

.2126 

.0164 

55.6 

.00606 

11.4235 

38 

9 ^ 

.17  35 

. 0575 

43.6 

.Q3''7S 

3 .7210 

39 

92 

.2402 

- 0953 

4 5.  8 

.057  30 

7.  7047 

4 0 

89 

.4562 

,1567 

‘f  4.  7 

.097o5 

7.0823 

41 

oti 

.1968 

. 0273 

52.  4 

.014  42 

9.8Q40 

42 

90 

.5890 

. 2224 

6 7.  6 

. 08  0 2 3 

18 .524? 

4 3 

91 

.439-J 

.1469 

46.  2 

.08790 

7.7143 

44 

QQ 

.0224 

. h163 

3 8.  2 

.28611 

5.4992 

45 

91 

. 59p5 

.1310 

43.  1 

.114  46 

6.8103 

4 fa 

93 

. 1911 

. 009^ 

6 o . 2 

. 00  3 72 

18.0922 

47 

92 

. Ob  37 

. GC22 

34.  0 

.00161 

4.7029 

4o 

92 

. 1911 

.0162 

5 9.0 

. 007  39 

13.0961 

49 

9 2 

. 1 796 

. 0 08  4 

4 4.  6 

. 005  19 

7.3425 

5C 

82 

. 03  72 

. 0 6 1 0 

6 0.4 

.00070 

12.1920 

51 

07 

. 0393 

. 0026 

6 3.  7 

.1)0105 

15.0222 

52 

88 

.2257 

. 0 39  0 

72.5 

.01197 

22.8657 

53 

ee 

.1^72 

, 01?H 

3 2.  •3 

. 0 09  3 0 

4.2557 

54 

87 

.4887 

.0440 

65.  6 

.01689 

16.34b6 

55 

08 

.114^ 

. OOSo 

4 3.  2 

. 005  44 

6.6033 

5b 

c8 

1.2928 

.4-^11 

4 7.  2 

.25397 

7.7429 

57 

87 

.2910 

. 1 09  3 

30.  3 

. 085  0 5 ■. 

3.7241 

58 

87 

. 212«< 

.0419 

2 2.  8 

.03511 

2.5327 

59 

87 

.9720 

. 3135 

2 0.  4 

.27857 

2.1970 

60 

8F 

. 39o0 

. 00o9 

2K„  0 

. 007  18 

3.2774 
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TABLE  A-7:  PRECIPITATION  ANP  CUBIC  RESPONSr.  SU.U-ACE  STATISTICS, 

TARGET  AREA  (CONTI NUEH) 


Band  No. 

No.  of 
Stations 

Mean 

,n  Precip.  y 

Precio. 

Var. 

100R“ 

MSi; 

61 

91 

. u 7 8 1 

. 0 068 

67.7 

.00317 

13.3733 

62 

91 

.h371 

.1086 

66.  6 

.04026 

17.9673 

63 

91 

.2502 

. 0 4l?  1 

6 0.6 

.01758 

13.8  214 

64 

80 

. 18  07 

. 0 754 

64.  8 

. 0299  7 

14.2912 

65 

81 

. 17C2 

. 0 26  3 

5 3.6 

.01228 

11.1733 

66 

94 

.0161 

. 0 00  9 

7 b.  5 

. 00-J  24 

30.4566 

67 

67 

. 0272 

.0010 

1 3,8 

. 000  96 

1.3667 

68 

86 

.2014 

. 0449 

56.  3 

.02135 

10.3536 

69 

o7 

. 1 J2.3 

. 0 135 

27.  4 

.010  91 

3.2  219 

70 

91 

.ul37 

. 0J04 

2 5.  3 

.000  36 

J.1336 

7 1 

5C 

, 0932 

. 0046 

48.  1 

.00  275 

8,2315 

72 

90 

. 0622 

• 0 u 2 7 

65.8 

.00104 

17.1197 

73 

91 

.06  29 

. 0 u 24 

6 3.  1 

.Cf’l  0 0 

15.4035 

7u 

91 

.1170 

. DO*.  1 

38.9 

.0028  0 

5.7201 

75 

96 

. 54  73 

, 0 335 

55,5 

.0  28  3 2 

11.9003 

76 

93 

.0394 

.0012 

5b,  8 

. ono  53 

12.1213 

77 

93 

.4737 

. 121  7 

51,8 

.3<'-5  06 

9.9042 

7 0 

93 

. 317^. 

. 0 448 

35.  0 

. 0 3324 

4,5472 

7fJ 

95 

. 3446 

. 0 778 

5 3.3 

.04016 

10.7997 

80 

9 6 

.119  3 

. 0295 

1 4.  8 

. 027  80 

1.664,3 

81 

96 

. 2573 

. 0 '24 

3 0.3 

.01722 

4 .lo31 

82 

94 

. C5  88 

. 0057 

21.  1 

. 004  98 

2.4945 

83 

94 

.12  0 5 

. 0297 

23.0 

. 025  13 

2.7916 

84 

9 3 

. 27  09 

• 0 j 0 6 

4 4.  7 

.01888 

7.4541) 

85 

94 

. 9 0 •<  0 

. iHQh 

57.5 

,089  7 1 

12.6234 

86 

86 

. 1006 

. UC59 

64.  2 

.00  2 34  ■ 

15.1757 

67 

90 

. 15  79 

. 0t33 

4 3.4 

.00238 

6,8233 

88 

90 

. 0301 

. 0 126 

26.  0 

.010  36 

3.1275 

89 

96 

. 1 1 316 

. 0074 

33.1 

. 0Q2v2 

4.  733  3 

90 

96. 

.2110 

. 02  36 

6 8.  0 

.008  36 

70.2723 
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TABLE  A-7:  PRECIPITATION  AND  ^UBIC  RESPONSE  SURFACE  STATISTICS 


TARGET  AREA  (CONTINUED) 


Band  No. 

No.  of 
Stations 

Mean 

,n  Precip.  y 

Precip. 
Var.  s2 

o 

o 

MSE 

u. 

<91 

94 

1 .0084 

. 3666 

41.  5 

,ZZ7Z7 

6.6270 

92 

. 6492 

.2032 

58.  7 

.09  3 04 

12.9739 

93 

91 

.2063 

. 0284 

45.  4 

.01723 

7.4734 

91 

.0170  ' 

. 0 u 06 

57, b 

.00029 

12.2033 

95 

94 

.1328  , 

. 0129 

5 9.  3 

.0058  0 

13.6173 

96 

9 3 

.6634 

' . 178  3 

75, b 

.048  25 

28.5423 

97 

93 

.4362 

.1397 

64.2 

. 055  36 

16.5645 

96 

93 

. 0940 

. 0077 

36.  4 

.00542 

5 .23  57 

99 

92 

. 8258 

. 1600 

45.  4 

.109 16 

7.56*45 

100 

69 

. 1082 

.0079 

« 

22.  4 

.00685 

2.5238 

10  1 

88 

i .0542 

. 002  0 

29.  8 

.00160 

3.6767 

102 

89 

. 0562 

.0017  ' 

49.  2 

.OUO  97 

8.499  3 

10  J 

95 

, .1838 

.0256  , 

4 2.  7 

.01632 

7.0‘»16 

lO** 

94 

. 07  83 

. 0077  ^ 

55.  4 

.00382 

11.6147 

ms 

89 

.4610 

.0693 

5 2.  ?■ 

.03689 

9.5975 

106 

85 

. 1624 

.0166 

65.  0 

.00650 

15.4504 

1 J7 

65 

. 2741 

. 0396 

7 7.7 

.00987 

29.0  76« 

I 


/ 
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TABLE  AS:  PRI.CIPITATICN  ANO  CUBIC  RESPONSE  SURPACE  STATISTICS, 

CONTROL  AREA 


Pand  No. 

No.  of 
Stations ,n 

Mean 

Precip.  y 

Precip. 

Var. 

lOOR^ 

MSF 

1 

IS 

.1013 

. C ICS 

7 4.  7 

. 007  7 1 

1 .6435 

2 

IS 

. Qn?6 

70.  2 

.OOS  47 

1. 31C5 

3 

16 

• j ^ u 0 

. u P 5 .3 

91.4 

.00071 

7.0945 

<• 

16 

. 0563 

. 0 01^ 

8fi.  6 

. 00  0 34 

5.2018 

5 

16 

. 1<.50 

. 0C3  4 

9 3.3 

.00141 

9.2296 

e 

16 

. 0975 

. 0039 

85.  9 

.00  1 36 

4.0593 

7 

16 

. 1125 

. 0‘‘8 

66.  0 

.00160 

4.0977 

8 

19 

.00  ~j6 

. 0 0 1 

5 7.5 

, 00  0 06 

1. 3520 

9 

19 

.010=; 

.0002 

6 9.  8 

.01011 

2.3136 

10 

19 

. j342 

.0005 

85.  1 

.000  15  . 

5.6911 

11 

19 

.18  0 0 

. 0029 

79.  7 

.0011" 

3.9241 

12 

19 

.2795 

.0114 

91.  9 

.001 85 

11.20  *)9 

1 3 

19 

.0734 

. 0045 

64.  5 

.00141 

5.4436 

1 4 

19 

.1000 

. 0 04  0 

55.  r 

.Of’357 

1.237’ 

15 

1 7 

.159*. 

. 0 073 

7 7.5 

.00377 

2.6553 

16 

17 

. 1229 

. 0 0’4 

6 2.  e 

. U"  288 

1. 313j 

17 

20 

. 892C 

. 0 360 

5 3.  6 

.0316  2 

1.2940 

16 

16 

.2383 

. 0 05  7 

75.  5 

.0''?95 

2 .7  321. 

19 

16 

. 546  1 

.0371 

78.6 

.01767 

2.  45'  2 

20 

If 

. 0.3  »5 

. 0u2  1 

64.  4 

.00169 

1.2074 

21 

17 

. 0647 

. U 1 w C* 

7 6.4 

. 005S6 

2.510.7 

22 

17 

. 1059 

. OCJo 

9 1.0 

.00057 

7 .8643 

23 

23 

. £365 

. 0 u4  2 

5 0.  S 

. or  3 5.7 

1 .47?: 

24 

23 

. 188  3 

. 0028 

3S.  7 

.00213 

1.8158 

25 

23 

. u6  7 0 

.0  015 

8 1.9 

. 0 0 v’  4 5 

6.5233 

26 

?3 

.0904 

. t J '♦  0 

o 4 . 1 

. 00243 

?.57j5 

27 

23 

. 10j9 

, 005  7 

91.  4 

.00182 

15.4134 

28 

23 

. 1570 

. 0066 

7 7.2 

.00255 

4.8916 

29 

23 

.M3'’ 

. 0 0 59 

75.  1 

.O'’?50 

*.  . 3528 

.1  n 

22 

.258.;’ 

.014? 

88.  1 

. Cr?97 

9. 3432 

Footnote;  Significance  probabilities  of  F-ratios  can  bo  approximated 

froTi  the  standard  F-table  with  9 and  (n-10)  degrees  of  freedom. 
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TABLE  A-8;  PRECIPITATION  AND  CUBIC  RESPONSE  SURFACE  STATISTICS 


CONTROL  AREA  (CONTINUED) 


Band  No. 

No.  of 
Stations 

Mean 

,n  Precip.  y 

Precip. 

Var. 

lOOR^ 

MSE 

31 

23 

. 1 1 iJ  0 

. 0046 

82.0 

.001 

6.5955 

6 c 

23 

. 0326 

. U003 

4 7.  7 

.00025 

1 .3152 

33 

23 

. 1217 

.0130 

72.  5 

.006  06 

3.8162 

34 

23 

. 2743 

. 0046 

70.  3 

.00239 

3.4250 

35 

23 

.1709 

.0023 

70.7 

.00136 

3.4933 

Zf- 

23 

.1083 

. 0022 

37.7 

.00231 

.8733 

37 

22 

. 2364 

. c062 

71.  0 

. 00  3 17 

3. 2605 

38 

22 

.2355 

. 0674 

92.  6 

. 008  73 

16.6791 

39 

^2 

. 4236 

. 0869 

65.6 

.0=;234 

2.5410 

40 

21 

.3338 

.0416 

07.9 

.00919 

8.3903 

41 

21 

. 1433 

.0108 

69.0 

.00608 

2.71'iS 

42 

21 

.3100 

, 0166 

76.9 

.00697 

4.0752 

43 

20 

.6280 

. 041  3 

6 3.  2 

. 028  86 

1.9122 

44 

19 

.4463 

. 027  2 

54.2 

.02492 

1.1823 

4 5 

19 

.1779 

.01^0 

7-..  8 

. 006  53 

2.9774 

4 6 

23 

.10  35 

. 0 092 

78.4 

.00335 

5.2523 

47 

24 

. 0317 

. 001  e 

6 8.  3 

. 00  0 34 

11  .7  033 

48 

Zi 

.14'^0 

. U025 

6 1.1 

.00  164 

2.2670 

49 

23 

. 20  17 

. 0125 

90.  7 

.00196 

14.1442 

50 

22 

.1495 

. 0C45 

86.3 

.001 09 

8.3912 

51 

20 

.0635 

. 0016 

84,  8 

. oon  -I? 

6.2201 

52 

20 

. 36  35 

. 0927 

7 5.  1 

.04382 

3.3559 

53 

21 

. 1952 

. 005  7 

94,  2 

.00060 

19.7731 

54 

22 

. 46ii0 

.017  0 

66.  6 

. 00993 

2.6623 

55 

22 

.1155 

. 0G47 

42.  0 

.00472 

.9660 

56 

22 

.7714 

. 1203 

46.4 

.11275  • 

1.1556 

57 

22 

.3414 

. 0621 

78.5 

.023  35 

4.8723 

58 

2? 

.1873 

. 0119 

90.  3 

.00203 

12.3431 

59 

22 

.6955 

. 0484 

62.9 

.03141 

2.2654 

60 

22 

. 0927 

. 0016 

83.  0 

.00  049 

6.515  0 
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TABI.I-;  A-8:  PRECIPITATION  .AND  CUSIC  RESPONSE  SORFACF.  STATISTICS, 

CONTROl.  AREA  (CONTINUED) 


Band  No. 

No.  of 
Stations 

^lean 

,n  Prccip.  y 

Precip. 

Var. 

lOOR^ 

MSF. 

f' 

6 1 

2 1 

. 1 70U 

. 0 145 

7 0.  7 

. 0G772 

2.9436 

62 

22 

.4745 

. 0733 

62.  2 

.043  78 

2. 1052 

b"' 

22 

. 4 -U  3 

. 0136 

5 4.  0 

.01496 

1 .5643 

64 

23 

. 4691 

. 0472 

7 3. 1 

.02150 

3.9192 

o5 

25 

. dyse 

. 029b 

7 7.7 

.01063 

5.315'v 

66 

25 

. 0372 

. 0 076 

37.  3 

.00154 

11,5002 

o7 

25 

. 0356 

.0013 

7 9.5 

.00042 

6.4707 

68 

25 

. 2o  20 

. 0 165 

61.  1 

.00498 

7.1477 

69 

24 

. 2096 

. 0U87 

5 3.-3 

. 00666 

1.7775 

7 0 

2 5 

. 0452 

. 0 009 

27.  7 

,00109 

.S39i 

71 

.2392 

. u 132 

61.7 

.003 12 

2.6841 

72 

24 

, 1536 

. 0026 

75.  4 

.001 00 

4 , 75o  ] 

73 

24 

.1400 

. un32 

60.1 

.00208 

2.3363 

7-* 

24 

.146:^ 

.0  112 

7 9.  1 

. 003  33 

5,399  3 

75 

25 

. 0744 

.0070 

8 7.7 

.00137 

11.9371 

/ b 

26 

. 0866 

. OO4I 

32.  -a 

.00426 

.SI  21 

77 

28 

. 28  54 

. 0 3 7 9 

64.  0 

,009  08 

10  .511  8 

74 

26 

.46/9 

. 0790 

6 4.6 

.018  78 

10,953b 

79 

2b 

. 42  4b 

.0153 

3 3.  6 

.01522 

1.021c 

80 

28 

. 19  t-4 

. r .'  . 3 3 

71.2 

.00791 

4.9429 

81 

26 

.10  39 

. onS2 

79.  0 

.00163 

7.5  4-'  i 

82 

26 

. 1304 

. j 120 

5 3,4 

.003  72 

2.0346 

83 

2b 

.1115 

.0130 

4 7.  4 

.01070 

1.502C 

84 

26 

. 1412 

.0464 

69.  1 

.00  3 09 

3.9  769 

85 

27 

. 3248 

. 0 49  7 

52.  3 

.03585 

2.>.t‘'8 

86 

24 

.2329 

. 0224 

9 3.0 

.00256 

20. 6-, 22 

6 7 

26 

. 1254 

. 0169 

84.2 

. 004  17  • 

9.4626 

88 

26 

. 1165 

. 0U5  1 

5 2.  8 

,00  3 '4  • 

1.9874 

60 

28 

. 0925 

.0110 

65.  1 

.006  79 

3 ,7  2o5 

90 

26 

, 0473 

. 0046 

74.6 

.00192 

5.2606 

-se- 


TABLE  A-8:  PRECIPITATION  AND  CUBIC  RESPONSE  SURFACE  STATISTICS, 

CONTROL  AREA  (CONTINUED) 


Band  No. 

No.  of 
. Stations 

Mean 

,n  Precip.  y 

Precip. 

Var. 

lOOR^ 

MSE 

91 

24 

. ?0?l 

. 0556 

6 0.7 

.03590 

2.4058 

92 

24 

. lliSO 

. 015  0 

74.  0 

.006  42 

4.4179 

93 

28 

.2636 

. 0415 

38.  2 

.030  43 

1 .2360 

9<* 

20 

. C593 

. 0034 

35.  5 

.00331 

1.1017 

95 

33 

.1133 

. 0 125 

6 0.6 

.00686 

3.9  3*»4 

96 

33 

. 30  24 

-.  0 55  3 

7 9.6 

.01572 

9.9459 

S7 

3 3 

. 26  Oo 

. 0409 

56.  8 

.02458 

3. 3644 

93 

3 3 

. 0652 

. 0010 

16.4 

.00187 

.5019 

99 

32 

.082“ 

. 0305 

25.  8 

.03190 

.8510 

lOU 

33 

. 10  39 

. 0 030 

81.2 

.000  79 

11 .0  09  6 

101 

33 

. 1252 

. 0 023 

5 3.  2 

.00150 

2.9058 

102 

34 

. 0694 

. 1)  0 1 2 

5 0.4 

. 000  05 

2.7086 

103 

?Q 

. W55 

. 0 16  0 

55.  6 

.010  44 

2.6425 

1 J4 

29 

.0714 

. 0097 

6 1.9 

. 005 

3.4274 

135 

31 

.50  19 

. 0331 

45.  3 

.025o2 

1.9335 

106 

25 

. 32  32 

. 0088 

69.  5 

.00428 

3.7919 

107 

24 

. 3438 

. 0153 

45.3 

. 01376 

1 . 2857 

-59- 


TA»LE  A-«:  ESTIKATbS  Of  COEFFICIENIS  Of  CUBIC  RECRESSICKS.  TARGET  AREA 
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TABLE  A-12:  MEANS  OF  PRECIPITATIONS  AND  TOEIR  VARIANCES  FOR  VARIOUS  GEOGRAPHICAL 
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